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^ ' Abstract 

The present paper deals with quantization of perfect fluid anisotropic cosmological 
models. Bianchi type V and IX models are discussed following Schutz's method of ex- 
pressing fluid velocities in terms of six potentials. The wave functions are found for several 
examples of equations of state. In one case a complete wave packet could be formed ana- 
qh lytically. The initial singularity of a zero proper volume can be avoided in this case, but 

^ ■ it is plagued by the usual problem of non-unitarity of anisotropic quantum cosmological 

models. It is seen that a particular operator ordering alleviates this problem. 

PACS numbers: 04.20.Cv., 04.20.Me 

1 Introduction 

The basic motivation behind a quantum cosmological model is two-fold. One is the fact that 
when the linear dimension of the universe reaches the Planck scale ( 10 _33 cm), the universe is 
indeed governed by a quantum picture. The second is the hope that a quantum description 
might be able to produce a singularity free birth of our universe. The quantum state of the 
universe is described by a wave function \1/ which is the solution of the Wheeler DeWitt equa- 
tion on a minisuperspace. For a comprehensive review, we refer to pQ. 

One critical problem in quantum cosmology is certainly that of a suitable choice of time 
against which the evolution of the universe is investigated. This is because the notion of time 
has different implications in General Relativity and Quantum Mechanics. Systematic attempts 
towards resolving this problem started in the 90 's [2]. The strategies also include a scenario 
where the notion of time plays no role whatsoever [3]. For a recent review of different strategies, 
see [I]. 



If matter is taken as a perfect fluid, the strategy adopted by Schutz[3in] becomes extremely 
useful as a set of canonical transformations leads to one conjugate momentum associated with 
the fluid giving a linear contribution to the Hamiltonian. The corresponding fluid variable thus 
qualifies to play the role of time in the relevant Schrodinger equation. In an ever expanding 
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model, the fluid density has a monotonic temporal behaviour, the time orientability is thus 
ensured. 

Schutz's formalism has been extensively used by Alvarenga et a/[7] for isotropic cosmological 
models. The investigation involves spatially flat, open and closed models. Matter content is 
taken as a perfect fluid with an equation of state p = a p. They found finite norm wave packet 
solutions of the Wheeler-DeWitt equations. One important finding is that some singularity 
free models could be constructed even without violating the energy conditions. A later work 
involves anisotropic Bianchi I cosmology [9j, which is the anisotropic generalization of a spatially 
flat isotropic model. Anisotropic models have problems, particularly that of non-unitarity. The 
Hamiltonian, although hermitian, is not self-adjoint. This is discussed in some detail in ref[9j. 

It is well known that a standard Copenhagen interpretation of a quantum cosmological 
model is not tenable. A "Bohm-de Broglie" interpretation normally does better. For a com- 
prehensive review, we refer to [$]. In the anisotropic case, the non unitarity leads to futher 
problems in the interpretation. As the norm of the wave function becomes time dependent, 
Bohmian trajectories are not conserved and the Bohm-de Broglie interpretation also becomes 
vulnerable [9]. It appears that the philosophy of interpretation might require a dramatic exten- 
sion so as to incorporate anisotropic cosmologies. 

In the present work Schutz's formalism is utilized to quantize Bianchi type V and type IX 
perfect fluid cosmological models which are the anisotropic generalzations of open and closed 
isotropic models respectively. In section 2 we discuss the formalism and use that to quantize 
Bianchi type V models. In section 3, a type IX model is quantized. Some concluding remarks 
are made in section 4. 

2 Schutz's formalism and quantization of Bianchi V cos- 
mological model 

The relevant action for gravity with a perfect fluid can be written as 

A= I d 4 x^f=gR + 2 [ d 3 xVhh ab K ab + [ d 4 x^g P (1) 

JM JdM JM 

where h a b is the induced metric over three dimensional spatial hypersurface which is the bound- 
ary DM of the four dimensional manifold M and K ab is the extrinsic curvature. Here units are 
so chosen that c = 16ttG — h is equal to one. P is the pressure due to the perfect fluid. The 
perfect fluid satisfies an equation of state 

P = ap, (2) 

where a < 1. This restriction stems from the consideration that sound waves cannot propagate 
faster than light. In Schutz's formalism [HE] the fluid's four velocity can be expressed in terms 
of six potentials. However, two of them are connected with rotation. In Bianchi V or IX models 
permit timelike geodesies which are hypersurface orthogonal, the rotation tensor vanishes, 
and one can write the four velocity in terms of only four independent potentials as 

u v =Ue tV + 9S, v ). (3) 
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Here h, S, e and 9 are the velocity potentials having their own evolution equations, where the 
potentials connected with vorticity are dropped. The four velocity is normalized as 



u v u v = 1. (4) 

Although the physical identification of velocity potentials are irrelevant for the formulation, h 
and S can be identified with the specific enthalpy and specific entropy respectively. This identi- 
fication facilitates the representation of fluid parameters in terms of thermodynamic quantities. 
The metric for the Bianchi V anisotropic model is written as 

ds 2 = n 2 (t)dt 2 - a 2 (t)dx 2 - e 2mx [ b 2 (t)dy 2 + c 2 (t)dz 2 } (5) 

where n(t) is called the lapse function. While a, b, c are functions of the cosmic time t, m is a 
constant. Bianchi type I model is recovered when m — 0. Eliminating the surface terms, the 
first and second terms of the action (1) give 



A g = I dt 



2 / ■ ■ 2 Joe 

— abc + acb + bca + 3n m — 
n\ a 



(6) 



where an overhead dot indicates a differentiation with respect to time t and A g is contribution 
of geometry to the action. So the gravitational Lagrangian density can be easily identified as 

2 • be 

C g = [abc + acb + bca) — 6nm 2 — . (7) 

If we now choose the three metric coefficients as 

a(t) = e A) , b(t) = e^^- , c(t) = e^+ + ^ (8) 



then we get 
and 



a(t) b(t) c(t) = e 3A) (9) 



C g = -2 [3 /9 - (/?'+ - $-) ] - 6eP°nm 2 . (10) 



e 

n 

Now Pq, (5 + and /3_ will be treated as the relevant variables in place of a, b and c. The 
corresponding conjugate momenta are then 

e 3/3o . 

Po = -12 A), 

n 

e 3/9 

P+ = 4 (/?+-/?_), 

n 

p 3/3o 

P- = -4 (11) 

The Hamiltonian of the gravity sector is now given by 

.2 ^2 



U g = nc~^ ( - |! + £± + 6e 4 *m 2 ^ , (12) 
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where momentum p_ is replaced in terms of p + using equation fill p . For the fluid part, the 
action can be written using thermodynamic relations for h and S [316]. The relevant equations 
are 

p = p (l + n) , h = l + U + P/p , TdS = <m + Pd(l/po) (13) 

where r, p, p and II are temperature, total mass energy density, rest mass density and specific 
internal energy respectively. Rewriting the third equation of ( fl3l) we get 

rdS = (l + n)d[ln(l + n)-alnpo] ■ (14) 

It then follows that, r = 1 + II and S = ln(l + II) — cdnpo- We can show that the equation of 
state takes the form 

P = %-—h^ a e-^ a . (15) 

(l + a) 1+1/a 

In a comoving system u v = (n, 0, 0, 0), and equation (J3J) yields 



A f = I dt 



n -^ e ^o + 9S) l+l l a e- s ' a 

(l + a 1+1/Q 



(16) 



Af being the fluid part of the action. As h > so (e + 6S) > 0. If we try the canonical methods 
used, for example, in [TO] the Hamiltonian for this action can be written in a very simple form 
with the canonical transformations 

T = p s e~ s p;^ , p T =p\ +a e s , e = e-(l + a )^ , p e = p e , (17) 

Pe 

along with p s = 6p e . Here p e = ps = and Cf, the Lagrangian density of the fluid, 
is the expression inside the square bracket of equation ( 1X6]) . The Hamiltonian for this perfect 
fluid can now be written as 

Uf = ne- Wo e 3{1 - a)t5o p T . (18) 

The advantage of using this method, i.e., using canonical transformations, is that we could 
find a set of variables where the system of equations is more tractable, while the Hamiltonian 
structure of the system remains intact. It also deserves mention that amongst the four velocty 
potentials mentioned, actually two are used, namely e and S. This is because e and h are 
related by 

u^e 4l = -h 

( see ref[5]) and one other, namely 6 can be settled using the normalization fl4]). 

The super Hamiltonian for the minisuperspace of this anisotropic quantum model is 



7x — H n -\~ 1~L 



9^ n f 

e -3/3 



n- 



24 



[-pl + Spl + lUe^m 2 + 2Ae 3{1 - a)fSo p T ] . (19) 



Here n acts as a Lagrange multiplier taking care of the classical constraint equation TL = 0. 
Using the usual quantization procedure [HI IT2] , we write the Schrodinger-Wheeler-DeWitt 
equation for our super Hamiltonian with the ansatz that the super Hamiltonian operator an- 
nihilates the wave function, 

-H |¥(A,,/3 + ,t)> = 0. (20) 
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There are attempts [12] to show that the classical Hamiltonian for a cosmological spacetime 
is zero, if one takes into account both the matter sector and the geometry sector together, 
like the present situation. But these attempts involves pseudotensorial calculations and the 
result depends on the minisuperspace chosen [14|. The problem perhaps stems from the fact 
that localization of energy in general relativity is not uniquely defined. Whatever be the status 
of the constraint ( )20|) in the most general case, we shall be using this following the standard 
practice. 

With p Xi — > —id Xi , pt —> idr, h = 1 equation (|20|) can now be written as 



tt(A),/3+,*) = 0. (21) 



In this equation T = t is the cosmic time co-ordinate if we choose the gauge n = e 3o ^° and this 
follows from Hamilton's classical equations as T = {T, Ti} = 3 "g . Now it must be mentioned 
that while constructing the Schrodinger-Wheeler-DeWitt equation (I2TI) we have considered a 
particular choice of factor ordering for po and e~ 3/3 °. We will discuss other choices of factor 
orderings and its consequences on our prime results in the subsequent section. We require that 
the super Hamiltonian must be hermitian, so the wave function ^ must satisfy these conditions 
[ISIIISIE]; 

^U^±oc = or ^'U^ioc = . (22) 

In order to solve for the wave function \l/ from equation fl2Tl) we employ a separation of variables 

as, 

v&(/3 ,/3 + ,T)=e- i£T £(/3 ,/3 + ), (23) 



we get 

H _ 3 |!| + 144 e 4A> m 2£ + 24E e 3(i-«)#,£ = o 

Op UP+ 



(24) 



where E is the separation constant. Further if we write 

£(/3o,/3+) = <KA)M/?+) , (25) 

we get 



1 90 _ 3 9_V_ 1U 4/3 2 , 24 E 3(l-a)A, = o 



(26) 



The solution for r] is 



v {P + ) = Ce^ + , (27) 



where C is the integration constant and k is again a constant of separation which has to be 
real so that the wave function is normalizable. The equation for <p is 

+ 144e 4/3o m 2 + 24Ee 3(1 - Q)/3o + fc 2 = O . (28) 
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For m = 0, this equation reduces to the corresponding equation for the Bianchi I model [9]. 
It is very difficult to solve this equation analytically for all allowed values of a. Here we will 
study the solution for some particular values of a. 
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2.1 a = — | (Distribution of Strings) 

With a = —1/3 equation ( 1251) becomes 

«9 2 ^ 



0fl 



+ 144e 4/3o m 2 + 24£e 4/3o + ^ 2 </> = O . (29) 



The solution of equation (129]) is known in terms of Bessel function and now we can write the 
wave function of equation (12"T]) as [T7] 



A; 

(A,, /?+, T) = e" ii?T [C X J& (V36m 2 + 6E e 2/3 °) + C 2 J^ (\/36m 2 + 6E e 2/3 ° ) ] , (30) 

where C\, C 2 are the integration constants. Now we can construct a regular wave packet 
superposing these solutions. In principle this can be easily done by considering the arbitrary 
integration constants to be suitable Gaussian functions of the parameters k and E ( see ref [9]). 
Defining q = y/36m 2 + 6E we can write the form of the wave packet as 

POO POO 2 t 

V wp = / f(k,q) q e i6m2T e~ l ^ T Jiu (ge 2 *) dk dq , (31) 

J k=—oo J q=0 

where f(k,q) = 2^ +1 e -M 2 -ik 2 an( j ^ an( j ^ are arD it rar y positive constants. The above 
integrals can be explicitly evaluated and the wave packet becomes 

= \/^77 e 9 e i6m T e e *J , (32) 



where = — In (A + ^). The norm of the wave packet is 

e 4 * * wp d(3 d(3 + = ^,/^eW , (33) 



6 

fOO POO 



6-7T 

2A 

oo ^ oo 



where w = arctan Clearly we can see that the norm is time dependent and hence the model 
is not unitary. This is not surprising as we know that the anisotropic quantum cosmological 
models are non-unitary [9]. The expectation value of any variable can now be calculated as 

lR N _ iz n ^ n P g; * WP dp dp + 

m n i: n P ^ P dp dp + ■ [6A) 

For (3i = (3o we find 

Wo) = ~l (In ^ + , (35) 



where B = \J X 2 + |^ and Z is a numerical factor ~ 0.58. Now 

e<*> = (a b c) 1 / 3 = (2Ae~')i (1 + Jg,)V* . (36) 

This describes the cosmological evolution of the spatial volume and this model predicts a bounce 
from a minimum volume universe with no singularity. Of course this result is facilitated by 
the choice of the gauge n = e 3a ^° as mentioned after equation(21). It is easy to see that the 
minimum volume is obtained at T = and the corresponding length scale is {2\e~ l )±. as 
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FiglirG 1 ! Plot of | \P | 2 against /3q and ft+ at a fixed time. The left panel shows the plot at T = and the right panel is the plot 
for T = 10. 

already mentioned, I = 0.58. So the linear size of the universe can be set in the model by 
choosing the vaule of A. For (3i = f3+, 

</?+> = -J (In 1 + /). (37) 

It deserves mention that (3+ is in fact a constant in time. The plot of the norm of the wave 

packet against (3$ and (3 + shows a well behaved pattern, and even at T = 0, does not show 

any sign of blowing up (see Figure [TJ). This is consistent with the fact that there exists a 

minimum of the proper volume and thus the pathology of a singular state of the universe can 

be avoided in this model. Furthermore, the plots are given for two different time, namely T = 

and T = 10 in some units, and the qualitative nature of the plots remain similar. 

Now let us have a closer look at equation fl33l) . As to = arctan ^ so we infer that the norm 

of the wave packet is time-dependent hence the model is non-unitary. We now investigate the 

i — 

time dependence of the norm (^ y^f as a function of time (T) in Figure [2]. It is observed 
that with the increase of T, the norm flattens, i.e., almost becomes a constant. This indicates 
that the problem of non-unitarity is somewhat diluted at least for a large time. It deserves 
mention that this feature may be facilitated by our choice of gauge n = e 3a/3 ° which makes T a 
constant. 

It seems important to investigate the same quantum cosmological model with factor order- 
ings other than that mentioned after equation (12 ip. The question of factor ordering arises as 
an artifact of the first term (— ^e _3/3 ° Pq) of the super Hamiltonian because in the quantization 
process we have operators assigned to each (3q and po- Earlier we have already studied the case 
where the ordering considered was e~ 3l3 °pQ. We now study two different cases of factor orderings 
where we consider the factor e~ 3/3 ° x as 

1. pi e" 3/3 ° 

2. p e" 3/3 ° po . 
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-40 -20 20 40 T 

Figure 2: Plot of Norm J^- e~) as a function of T as u> = arctan Here we have used A = .1 and 7 = 1. A suitable 
scaling is used to enlarge the figure. 



For the first case we found the solution for the wave function of the Schrodinger-Wheeler-DeWitt 
equation as 

k 

#i*(A),/9+,r) = e 3/3 ° e~ iET [C x (V36m 2 + 6E e 2/3 °) + C 2 J=ik (V36m 2 + 6E e Wo ) ] 

(38) 

and for the second case the same entity as 

iET J~kP+ 



^ 2 n d (/3 ,/3+,T) = e~ iET eVs p+ [C x (V36m 2 + 6E e 2/3 °) + C 2 (V36m 2 + 6E e 2/3 ° ) ] . 

(39) 

Here E' s and k' s are the separation constants for the respective partial differential equations 
concerned and C' s are the integration constants. If we compare the equations ( 1301) . ( 138|) and 
(|39|) we can clearly see that the only difference amongst them is marked by e T/3 ° where r = is 
for (|30|) . r = 3 for (|38|) and r = | for (|39|) . So now we consider 

._k_ R 

V(Po,P+,T) = e T/3 ° e" ii?T eW + [d (V36m 2 + 6£ e 2/3 °) + C 2 (v^m 2 + 6£ e 2/3 ° ) ] 

(40) 

as a solution of the Schrodinger-Wheeler-DeWitt equation where r is the footprint of different 
factor orderings considered. We re-write the super Hamiltonian for the minisuperspace of 
equation (TlQl) as 

T~i — 7~Lg + nt 

= _^ e -(3+ e +,)/3o po e e Po pQ g A) 

e -3/3o 

+ n-^- [3 p 2 + + 144 e 4/3 °m 2 + 24 e 3(1 - Q)/3 °p T ] . (41) 



8 



The relation between r, g and c; is given by 



and ( £?,?") can take values (0, 0), (0, -3) and (-3, 0). Now we re-examine the earlier results. The 
wave packet will also be modified and can be expressed as 



y wp = e TPo ^h e e e i6m T e e % . (43) 

An explicit calculation for the norm of the wave packet yields 



ny^Wpr / , rx /2\ 1+f A a , T 2 



2 , ,2 



e 4 * ^ dft = r (1 + - j J (^A 2 + - j eW . (44) 

We can also calculate the expectation values for (3q and j3 + but in that case we get the same 
results (equation ( 135]) and (1371)). So the expectation values of the metric functions will not 
be affected by the factor ordering as we would have the norm of the wave function in the 
denominator as well ( such as equation (34)). So the physical interpretation remain same. We 
now look back at equation (1441) . In Figure [3] we show the variation of the norm of the wave 
packet as a function of time for three different factor orderings. 

It is well known that the anisotropic models are plagued by the problem of non-unit arity. 
The problem appears to be generic. But our results show that a clever choice of operator or- 
dering, although unable to solve the problem, can provide a method for alleviating the problem 
to a certain extent. The operator ordering used in equation (30) ( and hence (35)) does better 
as figure (2) shows that the norm of the wave packet attains a constant value as T increases. 
For the two other choices, the expectation value does not lead to a constant value even for high 
values of T (figure (3)). 



2.2 a = 1 (Stiff Matter) 

With a = 1 equation (|28|) becomes 

d 2 ^ 



^ + 144 m 2 (f) + (24 E + k 2 )<j> = . (45) 

The solution is obtained in terms of Bessel function and we can write 

</> = C, (6m e 2 ^ ) + C,J (6m e 2ft > ) , (46) 

2 2 

where C3 and C4 are the integration constants. The wave function can now be written as 

(ft, + , T) = e-* ET e* v¥+ [C 3 J iV ^ (6m e 2 * ) + C 4 J (6m e 2 ^ ) ] . (47) 
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i e 8 t as a function of T as u) = arctan 4t ■ Here we have 



studied for three different factor orderings r = (Red), r = § (Blue) and r = 3 (Green). We have used A = 1 and 7 = 1. A 
suitable scaling is used to enlarge the figure. 



2.3 a = | (Radiation) 

With a = 1/3 equation (1251) becomes 

5 2 



+ 144 e 4/3 ° m 2 + 24 E e 2p °(b + k' z d) = . 



2/3o^, 1.2, 



(46 



The solution can be written in terms of Confluent Hypergeometric function (J 7 ) and Associated 
Laguerre function (£), 



__ e - ieme^o eim ^ T ( i(WH [ x + ^ ,i2me 2/3 ^ 



+ C 6 £ P (g ~ 2 t +fcm) ,^, ^12me 2/3 °) ] , (49) 
where C5, Cg are the integration constants. The expression for the wave function can be written 



as 



tt(A), T) =e~ iET e^ P+ e~ i&m ^ jk? ^ T x + ^ a2me 2/3 °) 

"' (fi - m+fcm) ^,2l2me 2 ' 30 ) ] . (50) 



3 Quantization of Bianchi IX cosmological model 

The Bianchi IX metric is written as 

ds 2 = dt 2 - a 2 (t)dr 2 - b 2 (t)d6 2 - [ b 2 (t) sin 2 9 + a 2 (t) cos 2 9 }d(f) 2 + 2a 2 cos 9 drd<f> . (51) 
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The gravitational Lagrangian density is written as 



r 2 ^ 2 2(32 a " ^9 



where (3 = a b. The Hamiltonian for this gravitational Lagrangian density becomes 



Ob 2 ^2 ^ 



n= w- p °-8 p * + w-- 2a ' (53) 

where p a and pp are the canonical conjugate momenta for the variables a and 0. Applying 
Schutz's mechanism [51 |6] and the canonical methods described in [10] we can evaluate the 
Hamiltonian for the fluid as 

U f = a a (3- 2a p T . (54) 
The super Hamiltonian for the minisuperspace of this model is 

a 3 2 a o cl 5 



Pa-*Pfi + ^- 2a + a "P PT ■ (55) 



8/3 2 ^ a 8 ^ 2f3 2 
We shall calculate the wave function of the Bianchi IX universe with a = 1. 



3.1 a = 1 (Stiff Matter) 

Using the quantization procedure as described in [HI [12] we write the Schrodinger-Wheeler- 
deWitt equation as 

_ ££i|AZ2 + ft T) - 2^(o, A T) + i^^Il = . (56) 

Now we have to apply the method of separation of variables for the solution of this equation. 
We write the wave function as 



which leads to 

With the separation 
we get 

and 



*(a,/3,T) = e-^£(a,/?) , (57) 
£(a,/3)=«(a)v(/9) , (59) 



^ <^ - ~ ^ = ( 60 ) 



a 2 — - 4a 4 u - 8(£ + fc)w = . (61) 
oa 2 
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Here k is the separation parameter. The solutions of (16"Uj) and (IBTj) are known in terms of 
modified Bessel functions of first (I) and second (/C) kind and can be written as 

v{/3) = (5 l/2 [ C 7 X v (24/3) + C 8 K v (*4/3) ] (62) 

and 



u[a) 



i l/2 [ C 9 T v , (ia 2 ) + C w Ky, (ia 2 ) ] , (63) 



where v = + 32A; and v' = + 32(i? + fc). Cj's are the arbitrary integration constants. 
So the final expression of the wave function is 

V{a,P,T) =e~ lET ^a~f3 [ C 7 X u (i4/3) 

+ C 8 /Q, (i4/3) ] [ C 9 2^ (m 2 ) + C 1Q K v , (ia 2 ) } (64) 



4 Discussion 

Bianchi type V and type IX perfect fluid models are quantized following Schutz's formalism. 
As the Wheeler-DeWitt equation for a general equation of state p = ap is too involved, some 
specific examples are taken up. In Bianchi type V models, the solution for the wave function 
for the universe comes out as a combination of Bessel functions for the case a — — |. But it 
deserves mention that the problem of non-unitarity exists in this model, the norm of the wave 
packet is indeed time dependent. A clever operator ordering can only alleviate the problem in 
the sense that the norm becomes a constant for a large time. But the problem is only alleviated 
and not eradicated! In addition to the operator ordering, the choice of a gauge as n = e 3a/3 ° 
also may have its say. 

The proper volume shows a bounce from a finite minimum and thus the singularity of a zero 
proper volume can be avoided. By a choice of parameter, the minimum length scale obtained 
can be tuned to a desired value, say the Planck length. However, this is obtained as a bonus, 
as the chief motivation of the work had been to check the merits and problems of the particular 
method of quantization for the anisotropic cosmological models. 

Only one case, namely a = — | has been solved here as an example. This is simply because 
this case could be integrated anayltically. For some other cases also the wave functions are 
obtained analytically. But the complete analysis could not be done as the wave packet could 
not be integrated. So the work is nowhere near a general investigation of the problem. However, 
so long there is not a complete resolution of the Quantum Cosmology, one has to look for various 
examples. 

So although potentially Schutz formalism gives us a beautiful way of handling the problem 
of the choice of time, it also leads to other problems, namely that of non unitatiry. This could 
be a generic feature of anisotropic models as indicated in [9]. The philosophy of interpretation 
also seems to be awaiting a new direction. 
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